Implicitly Restarted Arnoldi Methods and Eigenvalues of the Discretized Navier Stokes Equations

نویسندگان

  • R. B. Lehoucq
  • J. A. Scott
چکیده

Implicitly restarted Arnoldi methods and eigenvalues of the discretized Navier Stokes equations. Abstract We are concerned with nding a few eigenvalues of the large sparse nonsymmetric generalized eigenvalue problem Ax = Bx that arises in stability studies of incompressible uid ow. The matrices have a block structure that is typical of mixed nite-element discretizations for such problems. We examine the use of shift-invert and Cayley transformations in conjunction with the implicitly restarted Arnoldi method along with using a semi-inner product induced by B and puriication techniques. Numerical results are presented for some model problems arising from the ENTWIFE nite-element package. Our conclusion is that, with careful implementation, implicitly restarted Arnoldi methods are reliable for linear stability analysis.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Implicitly restarted Arnoldi with puri cation for the shift-invert transformation

The need to determine a few eigenvalues of a large sparse generalised eigenvalue problem Ax = Bx with semi-positive deenite B arises in many physical situations, for example, in a stability analysis of the discretised Navier-Stokes equation. A common technique is to apply Arnoldi's method to the shift-invert transformation, but this can suuer from numerical instabilities as is illustrated by a ...

متن کامل

Implicitly restarted Arnoldi with purification for the shift-invert transformation

The need to determine a few eigenvalues of a large sparse generalised eigenvalue problem Ax = λBx with positive semidefinite B arises in many physical situations, for example, in a stability analysis of the discretised Navier-Stokes equation. A common technique is to apply Arnoldi’s method to the shift-invert transformation, but this can suffer from numerical instabilities as is illustrated by ...

متن کامل

Krylov Methods for the Incompressible Navier-Stokes Equations

Methods are presented for t ime evolution, steady-state solving and linear stability analysis for the incompressible Navier-Stokes equations at low to moderate Reynolds numbers. The methods use Krylov subspaces constructed by the Arnoldi process from actions of the explicit Navier-Stokes right-hand side and of its Jacobian, wi thout inversion of the viscous operator. Time evolut ion is performe...

متن کامل

Some new restart vectors for explicitly restarted Arnoldi method

The explicitly restarted Arnoldi method (ERAM) can be used to find some eigenvalues of large and sparse matrices. However, it has been shown that even this method may fail to converge. In this paper, we present two new methods to accelerate the convergence of ERAM algorithm. In these methods, we apply two strategies for the updated initial vector in each restart cycles. The implementation of th...

متن کامل

A Harmonic Restarted Arnoldi Algorithm for Calculating Eigenvalues and Determining Multiplicity

A restarted Arnoldi algorithm is given that computes eigenvalues and eigenvectors. It is related to implicitly restarted Arnoldi, but has a simpler restarting approach. Harmonic and regular RayleighRitz versions are possible. For multiple eigenvalues, an approach is proposed that first computes eigenvalues with the new harmonic restarted Arnoldi algorithm, then uses random restarts to determine...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1997